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ABSTRACT 



Aims. Micrometeoroids (cosmic dust with size between a few /im and ~1 mm) dominate the annual extraterrestrial 
mass flux to the Earth. We investigate the range of physical processes occurring when micrometeoroids traverse the 
atmosphere. We compute the time (and altitude) dependent mass loss, energy balance, and dynamics to identify which 
processes determine their survival for a range of entry conditions. 

Methods. We develop a general numerical model for the micrometeoroid-atmosphere interaction. The equations of 
motion, energy, and mass balance are simultaneously solved for different entry conditions (e.g. initial radii, incident 
speeds and angles) . Several different physical processes are taken into account in the equation of energy and in the mass 
balance, in order to understand their relative roles and evolution during the micrometeoroid-atmosphere interaction. In 
particular, to analyze the micrometeoroid thermal history we include in the energy balance: collisions with atmospheric 
particles, micrometeoroid radiation emission, evaporation, melting, sputtering and kinetic energy of the ablated mass. 
Results. Low entry velocities and grazing incidence angles favor micrometeoroid survival. Among those that survive, 
our model distinguishes (1) micrometeoroids who reach the melting temperature and for which melting is the most 
effective mass loss mechanism, and (2) micrometeoroids for which ablation due to evaporation causes most of the the 
mass loss. Melting is the most effective cooling mechanism. Sputtering-induced mass loss is negligible. 

Key words. Meteors, Meteoroids - Atmospheric entry 



1. Introduction 

Micrometeoroids (hereafter referred to as /iMETs) with di- 
mensions between ~25 /jm and ~1 mm represent (by mass) 
the majority population of Solar s ystem minor bodi es at 
a distance of 1 AU from the Sun (jGriin et alJ 119851) and 
dominate the annual mass flux of extraterrestrial matter 
entering the Earth's atm osphere, i.e. 4 ± 2 x 10 7 kg/yr, 
(|Love fc Brownleel Il993tl . /iMETs can survive the inter- 
action with the atmosphere and reach the Earth sur- 
face, becoming micrometeorites. A large number of mi- 
crometeorites, > 10 5 , have been recovered in Greenland 



(iMaurette et al. _ 19871) and Antar c tica dMaurette et al 



19911 lEngrand fc Maurettd 119981: iDuprat et all l2007t 



Cordier et aLll201lUvan Gimieken et al.ll2012n " over the past 
three decades. Many of these recovered /iMETs are, how- 
ever, less altered than expected: unmelted micrometeorites 
account for between 10 and 30% for particles >100 /mi, 
~50% for sizes between 50 a nd 100 /im and up to 78% 
in the size range 25-50 /im (jGenee et al.l 120081 ). In the 
CONCORDIA collection, which comprises micrometeorites 
collected in 2000, 2002 and 2006 at Dome C, Antarctica, un- 



Send offprint requests to: 
como.briani@csnsm.in2p3.fr 



Briani, email 



gia- 



melted microme teorites represent between 33 and 36% of 
all the particles (jDobrica et al.ll2009L [2010l) . 

When interacting with the Earth's atmosphere, /iMETs 
are subjected to effects of different physical processes. 
These include collisions with atmospheric atoms and 
molecules, deceleration, temperature variations, ablation 
of /iMET mass and sputtering, and excitation and ion- 
ization of atmospheric atoms and molecules. Several mod- 
els have been proposed to analyze different aspects of the 
/iMET - atmosphere interaction. In many cases, models for 
the /zMET-atmosphere interaction focus on specific pro- 
cesses or are aim to explain particular meteor observa- 
tions. iFlvnnl (|l989h modeled the /iMET - atmosphere in- 
teraction using the U.S. Standard Atmosphere values for 
the atmospheric density instead o f the expo nential ap- 
proxim ation proposed by IWhipplel (|1950D and iFraundorl 
(|l980l ) , but assuming that t he /iMET mass remains con- 
stant during the interaction. lLove fc Brownleel (fl99l an- 
alyzed the heating experienced by /iMETs. Their model 
includes vaporization and ablation of the /iMET material 
but not sputtering, and they did not analyze how the differ- 
ent processes evolve during the /iM ET atmospheric flight. 
Campbel l-Brown fc Koschnvl ((2004) tuned their model to 
explain light curves of faint meteors (/iMETs in the size 
range 10 /im - 2 mm). In particular, they modeled the 



1 



G. Briani et al.: Simulations of micrometeoroid interactions with the Earth atmosphere 



/iMET fragmentation to explain sudden increases observed 
in the light curves and to ch aracterize the size of the grains 
that constitute the /iMETs. iRogers et al.1 (|2005f) compared 
the mass loss due to sputtering to that due to ablation. 
They assumed, however, that all the energy transferred to 
a /iMET by the collisions with atmospheric particles can 
be used by sputtering, rather than only a fraction of it, as 
should be the case because other physical properties (tem- 
perature incr ement, evaporation , etc.) depend on the col- 
lision energy. ISzasz et all (120081) added the treatment of 
sp uttering proposed by iTielens et ahl (|1994f) to the model 
of lLove fe Brownled (|1991l) to model /iMET head echoes 
rec orded by the UHF E uropean Incoherent Scatter facil- 
ity. IVondrak et al.l (|2008f) proposed a detailed model of the 
differential ablation of the /iMET elements. They included 
thermal ablation (i.e. vaporization and melting) and sput- 
tering and focused their analysis on the height at which the 
different /iMET chemical constituents (e.g. Na, Mg, Si, Ca, 
Fe) are released in the atmosphere. 

We approach the problem of the /iMET - atmosphere 
interaction in a holistic way, considering those physical pro- 
cesses that can be the most important and highlighting 
their relative contributions with different atmospheric en- 
try conditions and their evolution during the /iMET at- 
mospheric flight. This work is, however, not an exhaustive 
analysis of the /iMET-atmosphere interaction. Our goals 
are (a) to highlight the necessity of a general study (b) to 
propose a model that includes different physical processes 
and to analyze their respective contributions and (c) to de- 
lineate the most critical aspects of the problem. 



2. The model: properties of micrometeoroids and 
present-day atmosphere 

We treat /iMETs as homogeneous, constant density 
spheres, with radius t^met- Their initial radii are assumed 
to be between 25 /mi and 500 /im. The upper limit of 500 
/im is required by the assumption of isothermality: we as- 
sume that a /iMET is sufficiently small to achieve a uni- 
form internal temperature X^met through thermal conduc- 
tion. Actually, thermal gradients develop within /iMETs 
due to endothermic reactions [e.g. phase transitions, see 
iFlvnnl (Il995l)1 and are observed in recovered micromete- 
orites (iGengel 2006al) as w ell as in experimental simula- 
tions (|Toppani et al.1 [20011. But a treatment of /iMET 
internal temperature distribution is beyond the scope of 
the present work (a qualitative discussion is presented in 
section E]). Since we also assume constant emissivity, the 
/iMET size must be larger than the wavelength of the 
emitted and absorbed radiation. A /iMET generally ab- 
sorbs radiation in the visible and ultraviolet (VIS-UV) in 
the spctral range 200 - 800 nm and emits infrared (IR) 
radiation (at temperatures of ~1000 K) with wavelength 
?» 3 fiva (ICoulson fe Wickramasingell2003lh so the simula- 
tions terminate when r M MET < 5/im. Micrometeoritic struc- 
ture and compositions are particularly difficult issues. To 
treat the therm al processes occurring during their atmo- 
spheric passage (jGreshake et al.l fl998l : IGenge et al.ll2008lh 
we assume the /iMETs have the composition of carbona- 
ceous chondritic micrometeorites (which represent ~ 84% 
of the recovered micrometeorites, see e.g. lGengdl2006bl and 
iLevison et al.ll2009lh In particular, we assume that /zMETs 
are predominantly silicates. We choose the values of their 



thermodynamic parameters accordingly, as described be- 
low. As specified above, the /iMET density p^met is as- 
sumed to be constant during the atmospheric flight. We 
used two different values for p^met, 3 and 1 g/cm 3 , to 
span the range of density me asured for recovered microm- 
ete orites (IGenge et al. _ 20081) . S i milar values are adopted 



_ .Campbell-Brown fe Koschnvl (|2004[) and IRogers et al.l 
(|2005l) . 

The present day Earth atmosphere is represented by the 
altitude profiles of density and the conce ntration of ne utral 
components from the MSISE-90 model (|Hedirj Il99lh . We 
include six atmospheric species: N2, O, O2, Ar, He, H. We 
neglect the effects of charge accumulation in our models and 
therefore do not include atmospheric ions. Such processes 
are left for a future development. 

The atmospheric properties determine also how a 
/iMET moves in the atmosphere. Because the mean free 
path of an atmospheric particle, at heights between 60 and 
200 km (those of interest in this work), is always longer 
than the /iMET radius, we apply the free molecular flow 
regime. In other words, /iMETs always interact with sin- 
gle atmospheric atoms and molecules, and t here is not a 
gas layer that protects them from collisi ons (jSorasio et al.l 
[20011 : (Campbell-Brown fe Koschnvll200l . 



3. The model equations 

The two dimensional micromctcroid trajectory is described 
by the time dependent h (height from the Earth surface) 
and z (angle with respect to the zenith direction) . The equa- 
tions describing their evolution are: 



-77 - "^MET < 



dh 
dt 

and 



dz 
~dt 



smz / *Vmet 
1 9(h) 



"jiMET 



Ra 



(1) 



(2) 



where w m met is the /iMET velocity, g(h) is the gravita- 
tional acceleration at height h and is the radius of the 
Earth. The initial value, h — 400 km, applies to the present- 
day terrestrial atmosphere: the ambient density above this 
height is negligible. For the entry angle, different values 
were simulated (0°, 45° and 70°); for angles greater than 
70°, /iMETs tend to skip back to space. The initial value 
(vi n ) of the /iMET velocity v^met was left as a free param- 
eter, the observed ra nge is from 11.2 km/s to 72.8 km/s 
(jCeplecha et al.lll998D . 

Collisions with atmospheric atoms and molecules deter- 
mine the /iMET velocity v^met- We model this neglecting 
the thermal velocity of atmospheric particles (~370 m/s for 
T a tm = 160 K at h = 100 km) relative to i^met using the 
equation of motion: 



dv 



/iMET 



dt 



M, 



(iMET 



S> atm (M^MET + 9{h) COS z 



(3) 



where S = tit^met is the /iMET geometric cross section, 
M^uet is its mass, and p a tm(M is the altitude depen- 
dent atmospheric density. The drag coefficient V describes 
the momentum transfer efficiency in a single collision to a 
/iMET; its value is < T < 2 (r=2 for elastic collisions) 
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and we assume T = 1 (|Cam pbcll- Brown fe Koschnvll2004t 
iRoeers et al1l2005D . 

The core of our model is the description of how /iMETs 
acquire energy and of the changes thus produced. /iMETs 
acquire energy by collisions with atmospheric gas. Since 
vth <C h^met, the heating rate is given by: 



dE, 



coll 



dt 



A 

2' 



(4) 



where A, heat transfer coefficient, parameterizes the frac- 
tion of the kinetic energy of an incoming atmospheric par- 
ticle that is transferred to the target /iMET. Its value 
depends on the /iMET temperature an d on the /xMET 
Mach number A4 a MET = w^met / ^th (jHood fe Horanvl 
ll99lUMelosh fe Goldinll2008l) . However, in the limit where 
■Mumet ^S> 1, A ss 1 (see section [6] for a discussion) 
and we therefore assume A = 1, as do other recent stud- 



ies ( e.g. ICampbell-Brown fe Koschnvl 12004 IRoeers et al.l 
l2005h . Two other possible heating sources are s olar irra- 
diation ()Mosesl 119921 iMcAuliffe fe Christoull2006l ). and at- 
mosphere radiation. We did not include these sources in 
our energy balance because they are negligible for bodies 
as small as /iMETs. We estimated the contributions of so- 
lar and atmosphere radiation for r; n = 500 /im, v- ln — 11.2 
km/s, z- m = 45°, p^met = 3 g/cm 3 . The energy is 2 or- 
ders of magnitude larger than that gained from collisions 
at h = 262 km and 1 order of magnitude larger at h = 170 
km. However, this energy acquired by radiation does not 
provoke either a reduction in the size of the //MET size or 
a temperature increase. Energy acquir ed by radiat ion be- 
comes negligible, as in previous models (|Mosesl[r992D . when 
denser atmospheric layers are reached where the collisions 
dominate. In our example (rj n = 500 /^m, i>j n = 11.2 km/s, 
z- m = 45°, p^met = 3 g/cm 3 ) collisions start to dominate 
at h = 128 km and the maximum energy acquired by col- 
lisions overcomes that from radiation by 3 orders of mag- 
nitude. Thus, the pMET's temperature is regulated by the 
collisions: their peak temperatures and final sizes do not 
change when radiation effects are included. 

The rise in temperature due to collisional heating drives 
evaporation and melting of the /iMET, which reduce its 
mass and dimensions. Here the word "evaporation" refers 
to sublimation and that mass is assumed to be lost imme- 
diately. We assume that the evaporation begi ns as soon as 
mMETs start to acquire energy from collisions (iMoses 1992; 



ICampbell-Brown fe Koschny||2004l:rRogers et al.l2005l ) with 
the temperature dependent mass loss rate being mode led 
using the Knudsen-Langmuir formula (jBronshtenlll983l ): 



dt J 



= S 



27rfc B T' /i MET 



where /cb is the Boltzmann constant. Here the relation 
log 10 (p sa t) — A — B/T^MET (solution of the Clausius- 
Clapeyron equation) is used to evaluate the /iMET satu- 
rate d vapor pressure. We as sume A — 10.6 and B = 13500 
K f (|Love fe Brownleel 1 199 If )). We adopt /z = 45 amu for 
the mean atomic mass of the /iMET co nstituents that is 
appro priate for a chondritic composition (jLove fe Brownleel 

EmD). 

Melting also drives mass loss. At high temperatures, 
between ~1500 and ^2000 K, portions of the body un- 
dergo pass from a solid to a liquid and we assume that 
this melted mass is also completely lost. We note, however, 



that the very complex composition of chondritic /iMETs 
means that they have a range of melting points, so differ- 
ent parts within a //MET can begin melting while others re- 
main solid. We model melting using a term proportional to 
the he ating, as suggested by ICampbell-Brown fe Koschnvl 
(|2004f ). To enforce the condition that the rate of mass loss 
due to melting is different from zero only if T^met ~ Imeit 
= /iMET melting temperature, we employ a parametric co- 
efficient P S paii, defined as 



-Pspall = 0.5 



f spall = 0.5- 



arctan(X) 



(■MET 



If us) 



for T, 



7T 



^MET 



melt 



arctan[4(T M MET - T lus )} 



47T 



for T, 



/iMET 



melt ■ 



This coefficient was introduced by 

ICampbell-Brown fe Koschnvl ((2004), but in their work its 
value rapidly reaches 1 for T^met > 2melt> (i.e. all the 
heat is used to melt the /iMET mass). Since melting is 
only one of the processes included in the energy balance 
(@J), we adopt a form of P S paii that limits its value to 0.625 
(the two different expressions make P S paii a continuous 
function of T^met, to avoid problems in the numerical 
integration). The assumed value for the /iMET melting 
temperature is T mc i t = 1623 K. This is similar to values 
adopted in previous s tudies dBonnv fe Balageasl 119901 : 
lLove fc Brownleel 119911: iFlvnnl Il992t IScarsil 12004ft and is 
the melting point of silicates that are the dominant phases 
in /xMETs. The rate of mass loss caused by melting is thus 
given by 



dMfj,MET \ 

dt /, 



-Pspall / dE co \\ 
ffmelt V dt 



(6) 



where P^ mo it = 2.65 x 10 5 J/kg is the latent h eat of melting 
for stony meteoroids (|Love fe Brownleelll99lD . 

Sputtering is the ejection of atoms or molecules (tar- 
gets) by collisions with incident particles (projectiles) of 
sufficient kinetic energy. This mech anism has been prop osed 
to explain high altitude meteors (jBrosch et alJ 12001) and 
previ ous studies suggest that it can cause significant mass 
loss (jRogers et all 120051 ). We require the sputtering yield 
y, i.e. the number of target atoms or molecules removed 
by each projectile imp act. To calculate Y we ch ose a semi- 
empirical formula from lDraine fe Saloeterl ((19791) . again us- 
ing /Li = 45 amu as mass of the target particles and Uo = 
5.7 eV as their binding en ergy (this value is that in dicated 
for a silicate compound in lDraine fe Salpeterlll979l ). Using 
s for the total number of atmospheric species considered, 
the rate of mass loss by sputtering is: 



dM, 



/iMET 



(It 



sputt 



M-S^met riiYi 



(7) 



where rii is the number density of the i atmospheric 
species. 

In summary, evaporation, melting, and sputtering are 
considered in the ene rgy balance. As emphasized by 
lLove fe Brownleel (|l99lh . it is essential to include all three 
to properly evaluate the /iMET temperature. An important 
feature of our model is that the contributions of melting and 
evaporation are evaluated individually. The problem is to 
estimate the fraction 77 of the incoming energ y that is used 
for sputtering and not for the other processes. Pooova et al.l 
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(2007) proposed that particles removed by sputtering have 
a kinetic energy varying from 0% of the energy supplied by 
collisions (when the /iMET entry velocity is < 20 km/s) 
to 20% (when /iMETs enter the atmosphere at 70 km/s), 
and that the rest of the energy contributes to the rise of 
/iMET temperature that causes the emission of radiation 
and ablation. Based on this, and considering that a part of 
the incoming energy must be used to remove these parti- 
cles, we assume rj = 0.15 for every value of the /iMET entry 
speed. 

Beyond processes responsible for the mass loss, the rise 
of /iMET temperature causes intense emission of radiation. 
This cooling process is fundamental for the /iMET thermal 
balance, unlike larger meteoroids. As described in section 
[21 we assume /i METs have a constant emiss i vity, e = 0.9, 
as suggested in ICampbell-Brown fc Koschnvl (|2004ft . 

Finally we add a term f or the kinetic energy of the 
ablate d mass. As proposed in lCampbell-Brown fc Koschnvl 
pool , we assume that the ablated //MET mass leaves with 
a speed determined by the /iMET temperature T m met ■ The 
energy balance is thus given by: 

2 b P^m{ri)v piMET = c/w^MET — ^ h47rr MMET ec7i /iM ET+ 

( dMpMET \ „ ( dM^MET \ 
+-ncvap 77 + -"melt 77 + 

V dt y cvap v dt y molt 

. A 3 r-v. 1 / dMuMET\ SfceTuMET , aX 

+ 2^metL(w)+2 H (8) 

where a is the Stefan-Boltzmann constant, rji, m and rrii 
are the sputtering coefficient, number density and mass, re- 
spectively, of the atmospheric components. The expression 
for the energy used in sputtering takes into account that, 
depending on the collision velocity, sputtering is caused by 
only some of the atmospheric species (i.e. if Y is zero for 
the i th species, then r/i = 0). For the ther modynamic prop - 
erti es of /iMETs we adopt va lues used in lBronshtenl fl983) 
and lLove fc Brownled (|l99lft for stony meteoroids: c = 10 3 
J/ (kg K) is the /iMET specific heat, ff cvap = 6.05 x 10 6 
J/kg is the latent heat of vaporization. The variation of 
each term in equation [5] is followed throughout the numer- 
ical integrations so it is possible to see which are most re- 
sponsible for determining the /iMET fate and valuate how 
their roles change with different entry conditions and time 
during flight. 

4. Numerical integration and validation of our code 

Numerical integrations were performed using the RK4 
method with an adaptive step size. For each quantity (i.e., 
/iMET height, angle with respect to the vertical, speed, 
radius, mass and temperature) we required a minimum 
relative variation of 5 x 10 ~ 5 , and the maximum relative 
variation allowed was 10 -3 . If the relative variations were 
too small, the integration step size is doubled with an im- 
posed upper limit of 1 s; if the variations are too large then 
the calculated values are rejected and the integration step 
size halved with a lower limit of 10~ 10 s. We verified that 
the code reproduced the results of previo us studies. Several 
temperature vs time profiles are shown in lLove fc Brownled 



(11991ft . They use different initial temperatures and atmo- 
spheric properties than we do but we choose the difference 
between /iMET peak temperatures as parameter for a com- 
parison (i.e. the difference between the peak temperatures 
of /iMETs that have equal entry speed, angle and density, 
but with different initial si z e). Int roducing the set of equa- 
tions of lLove fc Brownled (|1991ft in our code, we obtain 
peak temperature differences that differ from the original 
ones for less than 9 %. A second comparison was made with 
iRogers et al.l(|2005h . who focused on the sputtering-induced 
mass loss. We computed the fractional mass lost by sputter- 
ing for /iMETs with p^met = 3 g/cm 3 and for /iMETs with 
P/iMET = 1 g/cm 3 , considering value of entry speed ranging 
between 20 and 70 km/s using the same set of equations of 
IRogers et all ()2005ft . The maximum difference between our 
results and those of IRogers et al.l (|2005h is 6%. 



5. Results 

Our simulations highlight how decisive entry conditions are 
B Tab. OQ summarizes the maximum values of V[ n that per- 
mits the survival of /iMETs with different initial radii, in- 
cidence angles, and densities. It is evident that 1) smaller 
/iMETs can have greater values of entry speed than greater 
/iMETs; 2) less dense /iMETs can sustain greater entry 
speeds than denser /iMETs and 3) /iMETs with large entry 
angle (i.e. far from the vertical) can survive higher entry 
speed than /iMETs with small entry angle (i.e. close to the 
vertical) . 



Table 1. Maximum values of entry velocity for which 
/iMETs are able to survive and reach the Earth surface. 





max Vi n 


(km/s) 




H 


P m met = 3 g/cm 3 


PpMET 


= 1 g/cm 3 
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= 70° Z in 


= 45° 
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21 


100 
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20 


20 


200 
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18 


18 


300 


11.2 12 


16 


16 


400 


12 


15 


15 


500 


11.2 


14 


14 



Figure [T] (all figures refer to /iMETs with p m met = 3 
g/cm 3 ) shows the radius evolution of a 50 /_tm /iMET for 
different values of v- m . In the case of z; n = 45°, for V[ n = 11.2 
km/s the final /iMET size is 28 |xm, i.e. the /iMET lost 
^82% of its initial mass, while for V[ n = 16 km/s the fi- 
nal /iMET radius is 7 /im (the mass lost is 99%). Another 
important feature evident in Fig. [1] is that /iMETs experi- 
ence most of their alteration within a limited height interval 



We pr esented preliminar y results of our model in a previ- 
ous work (Briani et al. 2007). In the present paper we correct a 
few errors in the model equations. We changed the expression 
for Papaii so that its maximum value is 0.65; we changed the 
equation for the rate of mass loss due to sputtering so that 
we now include sputtering caused only by those atmospheric 
molecules that collide w ith the /iMET with sufficient energy 
(jDraine fc Salpctcr 1979). In addition, we present here a more 
detailed analysis of how the considered physical processes evolve 
during the /iMET flight. 
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Fig. 1. Evolution of //MET radius with height assuming r u 
= 50 /im, z in = 45° and /9 m met = 3 g/cm 3 . 



(~15 km for the cases displayed in Fig. [IJ, corresponding 
to a duration of only few seconds. 

The same is true for the time evolution of //MET's veloc- 
ity. After a first phase in which the velocity remains almost 
constant, /iMETs are suddenly decelerated when they en- 
counter sufficiently dense atmospheric layers (Fig. [2]). The 
height at which /iMETs start to be decelerated is higher for 
larger entry radii and higher initial velocities but the time 
interval in which //METs change from their initial to their 
final velocity is generally short. A //MET with r m — 50 /mi, 
Ui„ = 11.2 km/s, z- ln = 45° and p m met = 3 g/cm 3 is slowed 
down to 1 km/s in ~4s, and then reaches its final velocity 
of ^100 m/s within ^30 s. In contrast, for rj n = 300 /mi, 
w in = 12 km/s, Zj n = 45° and p^met = 3 g/cm 3 , a //MET 
is slowed down to 1 km/s in only ~2.5 s. Those //METs 
with grazing trajectories or lower density experience lesser 
variation. For a very oblique entry, they decelerate more 
gradually. In particular, the //MET's velocity (and there- 
fore the rate of collisions) is reduced before the beginning 
of intensive ablation. Hence the maximum temperature is 
lower and the final radius is greater than for more nearly 
vertical entry. Similarly, //METs with low density values (1 
g/m 3 ) are decelerated at higher altitude where the atmo- 
spheric density is low. 

Typical thermal histories are shown in Fig. [3]A The 
almost impulsive temperature rise is followed by cooling 
with the whole process lasting only a few seconds. The ini- 
tial slight temperature decrease is due to efficient radia- 
tive emission and to the scarcity of collisions in the highest 
strata. //MET temperature would not decrease in the first 
phases of atmospheric entry if we included solar and at- 
mospheric radiative heating but, as previously discussed, 
this contribution has negligible effects on //MET peak tem- 
peratures and final sizes (section [3]) . At h ^100 - 120 km, 
the //MET temperature rapidly rises because of collisions, 
then after the maximum, T^met decreases very rapidly be- 
cause the collisional heating is reduced as //METs slow 
down and shrink. Taking the full width at half the max- 
imum (FWHM) temperature to compare the duration of 
the heating pulse, //METs in Fig. [3J^ suffer heating pulses 
of 5.78 s (for u m = 11.2 km/s), 5.25 s (for u m = 12 km/s) 
and 3.15 s (for w in = 16 km/s). Fig. [3]A_ shows that peak 
temperatures are reached earlier for higher entry velocities 
but the altitude of peak temperature is about the same: 
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Fig. 2. Evolution of pMET velocity with height (A) and 
time (B) assuming r ln — 50 /mi, Zi n = 45° and p m met = 3 
g/cm 3 . 



86 km for v- m = 11.2 km/s, 88 km for u; n = 16 km/s. Fig. 

and B show that higher entry speeds as well as larger 
initial radii imply higher peak temperature, but the dif- 
ferences are small. Comparing Fig. with Fig. [T] shows 
that the final radius is most sensitive to differences in entry 
speed. This is due to the melting process that both limits 
the maximum temperature values (Fig. [3J3) and generates 
significant mass loss. 

Fig. [3)3 shows temperature profiles for different val- 
ues of the initial radius (with v- m = 11.2 km/s). Again, 
higher peak temperatures are reached for larger ri n . Fig.[3j3 
shows the clear difference between two classes of profiles. 
For rin < 100 /mi the temperature curves have a regular, 
rounded shape. In contrast, for r- ln between 200 and 500 
/mi the top of the curves is due to melting of /iMET mate- 
rial. This is a general characteristic of our results. It is not 
linked only to different entry radii, but it appears for sev- 
eral combinations of entry conditions. These two classes of 
temperature profiles imply different final states for /iMETs. 
When the peak temperature is below the melting point and 
the temperature profile does not flatten, then the //MET 
survives after experiencing minor melting. These could cor- 
respond to the recover ed unmelted and pa rtially melted mi- 
crometeorites (see e.g. iGenee et al. 2008J). The other class 
of temperature profiles, with flatten peak near the melting 
point, correspond to /iMETs for which most of the mass 
is melted. The difference between these two classes can be 
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Fig. 3. Typical time variation of the /iMET temperature. 
A: results for r- ln = 50 /im, z- ln — 45° and p^met = 3 g/cm 3 . 
B: results for v in = 11.2 km/s, z- m = 45° and p^met = 3 
g/cm 3 . The truncated dashed and solid curves correspond 
to /iMETs that do not survive atmospheric passage. See 
text for discussion. 



seen also analyzing the history of the energy producing the 
temperature variation, i.e. the term cA^^met(^T^met/^) 
included in the energy balance. This term is displayed in 
Fig. |4j normalized to the incoming energy (equation [4]), for 
the same cases as Fig. [31 These curves cross zero around 40 
s, i.e. when the /iMET temperatures reach their peak value 
and then start to decrease. As shown in the insert in Fig. 
01 when /iMETs reach the melting temperature there is an 
interval during which no energy is used to alter their tem- 
perature. In contrast, curves corresponding to unmelted or 
partially melted /iMETs vary smoothly. 

Figure [5] shows the distribution of the energy that 
/iMETs gain from collisions in the different processes con- 
sidered in the energy balance. In the hrst part of the tra- 
jectory most of the energy goes into raising the /iMET's 
temperature (black, thin line). At high temperature, >1000 
K, ablative processes (evaporation and melting) and the 
kinetic energy of the ablated mass become important. 
However, the relative importance of the different processes 
change with different entry conditions. For the smallest 
/iMETs (Fig. the radiative losses dominate. This hap- 
pens in two phases. At the very beginning, not shown in 
Fig. \E\ when /iMETs suffer relatively few collisions, their 
temperature (Tj n = 330 K) is sufficient to produce signifi- 
cant emission. Then when the /iMET speed and radius are 




TIME (s) 

Fig. 4. Fraction of the absorbed energy used to change 
the /iMET temperature. These curves are proportional to 
dT^MET/dt so negative values mean that temperature is 
decreasing, positive values that it is increasing. The in- 
sert shows the transition from positive to negative values 
at around 40 s, coinciding with peak temperature and the 
start of /iMET cooling by intense ablation. 



greatly reduced, so is the energy acquired by collisions, the 
/iMET temperature is near its peak values, and the radia- 
tive loss dominates. 

Fig. [5K shows the typical case for a partially melted 
or unmelted /iMET. In the case shown, melting accounts 
for only about 0.1% of the incoming energy, while up to 
about 25% goes into evaporation. Sputtering is negligible, 
reaching about the same level as melting and terminating 
abruptly when /iMETs have slowed enough that the inci- 
dent kinetic energy of atmospheric particles is below the 
sputtering threshold. For the case shown in Fig. [5j3 abla- 
tion processes (evaporation and melting) are more impor- 
tant than in the previous case. There is a time interval 
during which both evaporation and melting consume more 
energy than radiative losses. Consequently, more mass is 
lost (see Fig. [5]) and the kinetic energy removed by the ab- 
lated mass is important for the energy balance. Sputtering 
is still unimportant because of the low initial velocity (v m 
= 11.2 km/s in Fig. [SJ3). Instead, for u; n = 14 km/s (with 
r; n = 500 /im, z- ln = 70° and /9 m met = 3 g/cm 3 ), we find an 
increase of about two order of magnitude, with the fraction 
of incident energy going into sputtering reaching ~10%. 
This is also found for /iMETs that do not survive atmo- 
spheric passage, as shown in Fig. [Sp; sputtering accounts 
for more than 10% in the energy balance, and melting and 
evaporation are more important than radiative losses. 

The analysis of the mass balance (Fig. H2) confirms what 
we find from the energy balance. In Fig. [6^ (unmelted or 
partially melted /iMET) melting accounts for ^20% of the 
total mass lost. For these /iMETs, mass loss is mainly from 
evaporation; sputtering is negligible. The situation is op- 
posite for the case shown in Fig. [5J3: melting is the main 
cause of mass loss, while evaporation driven mass loss is 
much lower (~16% for the case of Fig. [6)3) . In summary, 
sputtering is significant for greater entry velocities but, in 
our models, it is never the principal mechanism for mass 
loss. In the case of Fig. \EjP (relative to a /iMET completely 
destroyed) sputtering is responsible for ~0.1% of the total 
mass lost. Even in the most extreme case, i>; n = 70 km/s 
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Fig. 5. Variation of the energy going into different processes 
relative to the absorbed energy. A: unmelted or partially 
melted /iMET. B: almost completely melted //MET. C: 
completely destroyed /iMET. See text for discussion. Noise 
in the temperature variation lines is due to the numerical 
integration and occurs when the calculation is performed 
on very small step sizes (of the order of 10 -5 s, used when 
/iMET size, velocity and temperature vary significantly). 

and r in = 500 |xm, it accounts for a modest 1% of the mass 
loss. 



6. Discussion 

We have assumed constant values for both the drag T = 1 
(equation [3J and the heat transfer A (equation 3} coeffi- 
cients, a choice that is valid i n a fr ee molecular flow re gime 
[see e.g. ISorasio et all (|200lD and iPopova et all (|2007h ] at 
high Mmm- These coefficients are, however, dependent 
on the /iMET tempera t ure and the /iMET Ma ch number 
dHood fc Horanvilll991t livlelosh fc Goldfi] 12008ft . Thus, to 
test the sensitivity of our model to these parameters we 
performed simulations with different values of V and A for 
r in = 100 /im, = 12 km/s, z in = 45° and jo^met = 
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Fig. 6. Relative contribution to the mass loss of processes. 
The solid line (labelled "mass lost") shows the fraction of 
the total mass lost. A: evaporation dominates the mass loss, 
the final mass of this //MET is 47% of the initial mass. B: 
melting dominates the mass loss, the final mass is <1% of 
the initial mass. C: a /iMET that is completely destroyed 
by melting. 



3 g/cm 3 . Increasing T means the /iMETs are decelerated 
more effectively, so they are less heated by collisions, reach 
lower peak temperatures (peak temperature for T = 2 is 55 
K lower than for T = 1) and have larger final radii (r nn 
fi 7 /xm for r = 2, while r nn = 30 /mi for T = 1). In contrast, 
for r < 0.6 our test /iMET is destroyed; for V = 0.6 its peak 
temperature is 15 K higher than in the case T = 1 and the 
final radius is 6 /im. We find that /iMET peak tempera- 
tures and final radii are more dependent on A. Decreasing 
A from 1 to 0.1 yields lower peak temperatures, decreasing 
from 1604 K to 1255 K. The final radius is 99 /im for A = 
0.1, i.e. /iMET size is unaffect ed by the interaction with the 
atmos phere as also found by ICamobell-Brown fc Koschnvl 
(|2004f ). However, such low values of A are unlikely for the 
cases considered here. To estimate at what A the /iMETs 
experience the most intense alteration at peak tempera- 
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tures, we use the expression of iHood fc Horanvil (Il991h : 



A = 1 - 



7+1 



T, 



^MET/^air 



where 7=7/5 is the ratio of the specific heats of the terres- 
trial atmospheric gases and T a ; r is the atmospheric temper- 
ature. Since i?th is computed as Uth = V^te^air/^avg (&b 
= Boltzmann constant and m avg = mean mass of atmo- 
spheric particles) the value of A is determined by the ratio 
J/jMEt/^met- I n our simulations the maximum T^met is 
always reached before i^met is significantly reduced, and 
when w m met decreases, so does T m met- Thus, A remains 
close to 1. This is illustrated in Table 2. 



Table 2. Values of the heat transfer coefficient A for three 
test case /iMETs with z; n = 45° and p m met= 3 g/cm 3 . 
Simulations performed in the present-day atmosphere. 



Parameter 


Case 1 


Case 2 


Case 3 


r in (/im) 


50 


100 


500 


Win (km/s) 


16 


12 


11.2 


^TtMET (K) 


1606 


1603 


1622 


n P MET (km/s) at T^| T 


12.45 


9.65 


10.0 


ft (km) at T$&r 


88.0 


84.1 


79.5 


T air (h) (K) 


176.2 


181.4 


191.6 


i>th(ft) (m/s) 


392 


398 


409 


A 


0.973 


0.955 


0.958 



The coefficient 77, the fraction of the absorbed energy go- 
ing into sputtering, is also important but choosin g its value 
is ve ry difficult since previou s theoretical (e.g. ISigmundl 
Il98lh and experimental (e.g. iNakles! [2004) studies have 
treated sputtering as an isolated process. For instance, the 
mean mass of /iMET atoms and their binding energy are 
necessary parameters for specifying the sputtering yield Y, 
and its contribution to the energy balance. We note that 
a simplified treatment of the structure and composition of 
the /iMETs, as in this work, only approximates the sput- 
tering contribution and this process is probably more im- 
portant for less compact /iMETs than those treated here. 
Any treatment of the sputtering yield and heating hinges 
on several critical points: the evaluation of the different 
proces ses that cause th e removal of a particle from a solid 
body (|Mav et al.l l2000D . the spu ttering yield dep endence 
on the projectile incidence angle (jJurac et al.lll998lh which 
particles are ejected from th e target (atoms? molecules? 
both? See lTielens et al. 19941) and th e speed at which they 
leave the target (|Popova et al.ll2007t ). ICoulsonl (|2002h as- 
sumed a Maxwellian distribution for the speed of sputtered 
particles. This produces a momentum transfer to /xMETs 
about 11 times greater than the case of collisions alone, 
with consequent significant decreases of the /iMET speed, 
peak temperature and ablation. In our model, we do not 
include any impulsive reaction sustained by /iMETs from 
removed fragments since we assume isotropic ablative mass 
loss, although we do account for their kinetic energy in the 
energy balance. 

Finally the atmospheric component number density also 
affects sputtering. We find that sputtering-induced mass 



loss is always a small fraction (< 0.01) of the total am ount. 
This contrasts with the results of iRogers et al.l (|2005l ). who 
find that sputtering accounts for up to about 40% of the 
mass loss for /iMETs with p^met = 3 g/cm 3 and up to 
half for /iMETs with p m met = 1 g/cm 3 . We suggest that 
the difference betw een our results and those obtained by 
IRogers et al.l (|2005|) is due to their assumption that all the 
energy transferred to /iMETs by collisions with atmospheric 
particles produces sputtering, rather than only a fraction of 
it (i.e. they assume 77 = 1) and to their overestimate of the 
number densities of the atmospheric components that are 
systematically two order of magnitude greater than those 
we obtained from the same atmospheric model (the MSISE- 
90 model). 

Our model assumes that /iMETs are homogeneous 
and that all have the same composition. Recovered mi- 
crometeorites are, however, complex aggregates of different 
phases (e.g. silicates, sulfides, carbonaceous material) and 
show a wide range of compositions. More realistic /iMETs 
could be modeled assuming different proportions of the 
constituent phases. This would require modeling /iMETs 
of different densities and with different melting temper- 
atures. Our results show that a lower /iMET density fa- 
vors their survival (see Table [T]) because they are decel- 
erated at higher altitude, in less dense atmospheric lay- 
ers. Qualitatively, we expect that lower melting temper- 
atures correspond to those /iMETs that are mainly com- 
posed of more fragile phases, i.e. phases easily removed by 
the atmospheric frictional heating, as carbonaceous mate- 
rial and that higher melting temperatures correspond to 
those mainly composed of refractory phases, such as sili- 
cates, t hat are less affected b y the atmospheric frictional 
heating. IVondrak et al.l (|2008f ) modeled the differential ab- 
lation of the /iMET constituents, but without consider- 
ing heterogeneities and retaining the hypothesis of uni- 
form /iMET temperature. We too assumed that /iMETs 
have a uniform temperature at each in stant of their at- 
mosp h eric flight, as in previous mod e ls (|Lpve fc Brownled 
119911: ICampbell-Brown fe Kosclmvl [2004J). The prob- 
le m of internal t e mperature g radients was treated 
bv iSzvdlik fc Flvnnl (fl992h iFlvnnl (fl995h JSzvdlik fc Flvnnl 
(jl997ftlGenge fc Gradvl (|2000t ). However, in all these stud- 
ies, the /iMET radius is kept constant as are the thermal 
conductivity and heat capacity (i.e. they assume homo- 
geneous composition). A model for the simulation of the 
atmospheric passage of heterogeneous, multiphase /iMETs 
and of their internal temperature distribution would be a 
major advance. 

7. Conclusions 

To summarize, we list here our major conclusions. For any 
initial radius, a higher entry velocity and a more nearly 
normal incidence yield more severe heating for /iMETs, and 
therefore greater destruction rate. For fixed entry velocity 
and angle, the larger the radius, the greater the amount of 
mass lost. For all the /iMETs simulated in this work ra- 
diation emission is the most important process at greatest 
altitude, before /iMETs encounter atmospheric layers suf- 
ficiently dense to activate ablation processes (evaporation 
and melting). Melting is the most effective process that 
limits the increase of the /iMET temperature and produces 
significant mass loss. In all cases, the contribution of sput- 
tering is found to be negligible in both the energy and the 
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mass balance but could be more important for more fragile 
/iMETs than those simulated here. For the smallest /iMETs 
simulated in this work (r- m = 25/xm), radiation is the dom- 
inant cooling process during all their atmospheric flight, 
allowing the survival of these small /iMETs. Ablative pro- 
cesses are more important for larger /iMETs, i.e. these lose 
a larger fraction of their mass than small /iMETs and only 
low entry velocities and grazing entry angles insure survival 
during atmospheric passage. 

Our model shows that, among /iMETs that are not com- 
pletely destroyed, two types can be distinguished based 
on their temperature evolution. The first group is those 
/iMETs that show smooth temperature curves, i.e. that are 
not completely melted during atmospheric flight and there- 
fore reach the surface as unmelted or partially melted mi- 
crometeorites (most of the mass lost by these /iMETs is 
due to evaporation). The second comprises /iMETs with 
flattened temperature curves near the maximum tempera- 
ture. As we have shown, these differences are due to the 
importance of the melting process in the energy and the 
mass balance equations. 

We view this model as a first step in a broader study 
of /iMETs-atmosphere interactions. The necessary refine- 
ments require a broadly multidisciplinary approach that 
combines the physics with laboratory studies of the mate- 
rial properties from chemistry and mineralogy, for the alter- 
ation sustained by /iMETs in flight, with detailed numerical 
hydrodynamic simulations of the composite materials typ- 
ical of real samples. Finally, the analysis of the biochem- 
ical aspects, i.e. related to the organic molecules present 
in /iMETs, could improve the understanding of what role 
/iMETs played in the emergence of life on Earth (and, by 
extension, extrasolar planets). For example, although its is 
beyond the scope of the present paper, this type of modeling 
can be applied to the case of the interaction of /iMETs with 
the primordial terrestrial atmosphere to ass ess the role o f 
/iMETs as conveyors of pre-biotic molecules (Briani 2010). 
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